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ABSTRACT 

We investigate the connection between starburst and active galactic nucleus (AGN) activity by comparing 
the 3.3 /im polycyclic aromatic hydrocarbon (PAH) emission with AGN properties. Utilizing the slit-less 
spectroscopic capability of the AKARI space telescope, we observe moderate-luminosity Type I AGN at z^O.4 
to measure global starburst activity. The 3.3 /im PAH emissions are detected for 7 out of 26 target galaxies. 
We find no strong correlation between the 3.3 /im PAH emission and AGN luminosity in the limited range 
of the observed AGN luminosity, suggesting that global star formation may not be closely related to AGN 
activity. Combining our measurements with the previous 3.3 /im measurements of low redshift Type I AGN in 
the literature, we investigate the connection between nuclear starburst and AGN activity. In contrast to global 
star formation, the 3.3 /im PAH luminosity measured from the central part of galaxies correlates with AGN 
luminosity, implying that starburst and AGN activity are directly connected in the nuclear region. 

Subject headings: accretion, accretion disks - black hole physics - accretion - galaxies:active - galaxies: nuclei 



1. INTRODUCTION 

The correlations of black hole mass with galaxy proper- 
ties imply a close link between black hole growth and galaxy 
evolution i Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
iGultekin et al. 2009; Woo et al. 2010). Numerous theoretical 
and observational studies have been dedicated to revealing the 
nature of the connection (e.g., Kauffmann & Haehnelt 2000; 
Springel et al. 2005; Croton et al. 2006; Woo et al. 2006, 
2008; Hopkins et al. 2007; Merloni et al. 2010; Bennert et al. 
2011). 

Probing the direct connection between starburst and ac- 
tive galactic nuclei (AGN) can shed light on how black hole 
growth and bulge growth are inter-connected. However, it 
is challenging to observationally measure both the starburst 
rate and AGN activity for given galaxies. Optically-identified 
broad-line AGNs are the best objects for measuring black hole 
mass and AGN accretion rate. On the other hand the star for- 
mation rate of host galaxies is difficult to measure since AGN 
flux dominates in the UV continuum and hydrogen recombi- 
nation lines, i.e.. Ha, which are generally used as a star for- 
mation indicator in non-AGN galaxies. In the case of Type II 
AGN, it is relatively easy to measure star formation rate while 
black hole mass and accretion rate are difficult to determine 
owing to the dust obscuration of accretion disk and broad-line 
region. 

The polycyclic aromatic hydrocarbon (PAH) emission fea- 
tures emitted by small carbon grains are one of the star for- 
mation indicators. Thus, PAH emissions have been used to 
constrain star formation activity in the host galaxies of Type 
1 AGN (e.g.. Ogle et al. 2006; Schweitzer et al. 2006; Shi 
et al. 2007). Based on the ground-based and space obser- 
vations, various PAH emission features, i.e., 3.3, 6.2, 7.7, 
8.6, 11.3, and 12.7 nm lines, have been detected in many 
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local Seyfert 1 galaxies and distant QS Os althoug h not all 



AGN s how PAH emission features (e.g. , Clavel et al. 200C; 
Rodriguez- Ardila & Viegall l200l ifinanishi & Wada 2004 



Netzer et all 120071; iLutz et 31.1120071 120081; iShi et al] [20091; 
Imanishi et al.ll20ld;IOi et alj|2010l;llmanishi et 31^2011 1)! 



Among PAH-detected AGN, there seems to be a correlation 
between nuclear starburst and AGN activity. Using a sam- 
ple of low -redshift Palomar-Green (PG) QSOs observed with 
the Spitzer Infrared Spectrograph (IRS), Netzer et al. (2007]) 
showed that 7.7 /im PAH luminosity increases with AGN con- 
tinuum luminosity at 5100A (L5100), indicating the starburst- 
AGN connection. High-redshift QSOs seem to follow the 
same tren d, extending the co rrelation to much higher luminos- 
ity range (iLutz et al.l 120081) . The starburst-AGN connection 
is supported by a number of theoretical studies, demonstrat- 
ing that both nuclear starburst events and AGN activity can 
be triggered by gas inflow through galaxy interactions (e.g., 
Kauffmann & Haehnelt 2000; Robertson et al. 2006; Ciotti & 
Ostriker 2007; Hopkins et al. 2009). 

For lower luminosity AGNs it is less clear how AGN ac- 
tivity is connected to starburst since AGN power is more dif- 
ficult to measure owing to increasing host galaxy contribu- 
tion in lower lum inosity AGN. Based on t he ground-based L- 
band slit spectra. Umanishi & Wadal ( 12004 ) measured the 3.3 
/im PAH luminosity of local Seyfert 1 and 2 galaxies, showing 
a positive correlation between starburst and AGN activity al- 
though the correlation is not tight with a large scatter In their 
study, radio luminosity and N-band luminosity within 1 .5" of 
the host galaxies are used as an indicator of AGN power. Us- 
ing 2-10 kev X-ray luminosity, Watabe et al. (2008) showed a 
similar trend between AGN power and the 3.3 /im PAH lumi- 
nosity for a larger sample of local Seyfert 1 and 2 galaxies. 

Most previous observational studies on PAH features with 
ground based long-slit spectroscopy or Spitzer IRS observa- 
tions have probed the inner part of galaxies, indicating a con- 
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nection betw een nuclear starburs t within a few kpc and AGN 
activity (e.g. JWatabe et al.ll2008l) . In contrast, whether global 
star formation over entire galaxies is related with AGN activ- 
ity it is not well studied since most previous studies, partic- 
ularly for local and low-redshift galaxies, have used a rela- 
tively small aperture to extract PAH emission features, mea- 
suring nuclear PAH emissions. Based on the spatially re- 
solved mid-IR spectroscopy of nearby active galaxies, sev- 
eral studies have reported that PAH emissions are mainly from 
the disk (e.g., Laurent et al. 2000; Le Floc'h et alJi 2001). sug- 
gesting that the global PAH luminosity measured over entire 
galaxies may be much larger than nuclear PAH luminosity. 

In this paper, we investigate global and nuclear star for- 
mation in the host galaxies of Type 1 AGN, using the 3.3 
/xm PAH emission measurements based on our new obser- 
vations with the AKARI space telescope and collected data 
from the literature. The AKARI telescope provides a unique 
capability of slit-less spectroscopy to measure global star for- 
mation over entire galaxies (e.g., Imanishi et al. 2010). This 
paper is organized as follows. In ^ we describe sample se- 
lection, observations and data reduction of the AKARI obser- 
vations. We present 3.3 /im PAH emission measurements in 
^ Main results on the starburst-AGN connection are pre- 
sented in ^ Discussion and conclusions are followed in § 5. 
Throughout the paper, we assume a Hubble constant of H,, = 
70 km S-' Mpc-\ f^A = 0.7, and = 0.3. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample Selection 

We selected a sample of 27 moderate-luminosity AGN 
(L5100 0.1-3 X lO-*"* erg s'^) at 0.35 < z < 0.36 from the 
Sloan Digital Sky Survey Data Release 7 to study star forma- 
tion activity based on the PAH feature at 3.3 /im. The sam- 
ple was initially selected for investigating the evolution of the 
black hole mass-stellar velocity dispersion (Mbh -cr*) rela- 
tion. D eta ils of sample sele ction can be found in IWoo et akl 
( I2OO6I) and lTreu et alldlOOT i). In summary, we selected broad- 
line AGNs with the H/3 equivalent width (EW) > 5 A to 
properly estimate black hole masses based on the kinemat- 
ics of broad-line region gas. The specific redshift range was 
selected to avoid sky emission lines on the stellar absorption 
features, from which stellar velocity dispersions were directly 
measured. The spectral energy distribution of 24 objects in 
the s ample was in detail investigated by (Szathma rv et al.l 
I2OI lb based on the multi-wavelength data including Chan- 
dra X-ray observations, HST imaging, and Spitzer IRAC and 
MIPS imaging as well as the archival data of GALEX and 
Two Micron All-Sky Survey (2MASS). From the spectral en- 
ergy distribution (SED) fitting analysis, AGN SED was sep- 
arated from galaxy SED, leading to accurate measurements 
of AGN bolometric luminosities. Black hole mass of each 
AGN was estimated by the single-epoch virial method, which 
combines the line width of H/3 and the optical luminosity at 
5100A(L5ioo), utilizing the size-luminosity relation of the re- 
verberation sample (Bentz et al. 2009) and the m ost recent 
calibration of the virial factor by Woo et al] ( 1201 0) (see Ta- 
ble 1). The black hole mass of the sample covers an order 
of magnitude, i.e., 7.6 < log M^hIMq < 8.7. Combining 
black hole mass estimates with the bolometric luminosities 
measured from the SED fitting, we obtained the Eddington 
ratios, ranging from a few to ^50% of the Eddington limit. 
The basic properties of the sample are listed in Table [T] 



2.2. AKARI Observations 

To detect the PAH features at the rest- frame ^3.3 /im, 
we carried out spectroscopic observations using the In- 
frared Camera (IR C) on board the AKARI space telescope 
jOnakaet al.ll2007l) . The field of view of IRC is roughly 10' 
X 10' with a spatial scale of 1.45"/pixel. The full-width-at- 
half-maximum of the point spread function is 3.2 pixels, cor- 
responding to 4.64". The NIR grism mode (NG) was adopted 
for most of the observations, providin g spectral resolution of 
R ^ 120 at 3.6 /im for a point source jOnaka et aTll2007h . A 
unique capability of IRC is the slit-less spectroscopic mode. 
With r by r aperture used for pointing mode, the photons 
from the entire galaxies in our sample can be detected. Thus, 
the measured 3.3 fim PAH luminosities represent the global 
star formation. 

The exposure time for each AGN host galaxy was deter- 
mined in order to achieve signal-to-noise ratio S/N > 10 in 
the continuum. We assigned 4 pointings (6 minute expo- 
sure per pointing) for each target based on the exposure time 
calculation using the K-band magnitudes from the 2MASS 
archive. Owing to the complexity of scheduling of the sun- 
synchronous polar orbit of the AKARI, a couple of objects 
were observed with less than 4 points. In summary, we ob- 
tained a total of 93 pointings for 27 target galaxies (see[T]|. 

2.3. Data Reduction 

We carried out data reduction using the IRC spectroscopy 
pipelin^l The data reduction procedures with the IRC spec- 
trosco py pipeline is described in detail by lOhvama et alj 
(2007). In summary, the reduction steps include standard pro- 
cedures: dark subtraction, linearity correction, flat-fielding, 
background subtraction, extraction of 2-dimensional spectra, 
wavelength calibration, and flux calibration. Before extract- 
ing 1 -dimensional spectra, we apphed additi onal cosmic ray 
removal using L.A. Cosmic dvan Dokkumf 2001) to remove hot 
pixels, which significantly increased in Phase 3 (post liquid- 
He warm mission cooled by the onboard cryocooler). To ex- 
tract 1 -dimensional spectra, we us ed five pixels along the spa- 
tial direction as recommended by lOhvama et alj (120071) . As 
a consistency check, we compared the spectra extracted from 
larger aperture sizes and found no difference in the extracted 
spectra. Finally, we combine the extracted 1 -dimensional 
spectra from each pointing while applying additional sigma- 
clipping. In order to improve signal-to-noise ratio, we apply 
three pixels binning along the wavelength direction. We were 
not able to extract a spectrum for one object (SIO) due to the 
low quality of the observed spectral images. Thus, this target 
is excluded for the following analysis. 

3. 3.3 ^im PAH MEASUREMENTS 

3.1. Individual galaxies 

In Figure [T] we present the reduced spectra of 7 galaxies 
with the 3.3 fim PAH emission detection. Due to the low 
signal-to-noise ratio of the observed spectra, it is necessary to 
assume an intrinsic profile of the 3.3 /im PAH emission line 
in order to confirm the detection of the line and measure its 
strength. In measuring the line strength of the PAH emission 
lines, several different fitting schemes with various line pro- 
files, i.g., Gaussian, Lorentzian, Drude profil es, and Spline 
functions have been devised. For example, iGalliano et alj 
(.2008.) applied Spline functions and Lorentz profiles for fitting 

' |http://www.ir.isas.jaxa.jp/ASTRO-F/Observation/DataReduction/IRC/| 



3 



TABLE 1 

Properties of the AKARI sample: Type I AGN at z~0.4 



Name RA (J2000) DEC (J2000) redshift \oghtoi ALsioo.mir fnuc log Lx log Mbh pointings 

(erg/s) (erg/s) (erg/s) (Mq) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



SOI 


15 39 16.23 


+03 23 22.1 


0.3592 


45.23 


0.72 


0.29 


44.48 


8.20 


1 


S02 


16 11 11.67 


+51 31 31.1 


0.3544 


45.13 


0.34 


0.22 


44.60 


7.98 


4 


S03 


17 32 03.11 


+61 17 52.0 


0.3583 


45.50 


1.64 


0.39 


44.95 


8.28 


4 


S04 


21 02 11.51 


-06 46 45.0 


0.3578 


45.27 


1.33 


0.36 


44.55 


8.44 


4 


SOS 


21 04 51.85 


-07 12 09.4 


0.3530 


45.44 


1.85 


0.47 


44.84 


8.74 


4 


S06 


21 20 34.19 


-06 41 22.2 


0.3684 


45.05 


0.51 


0.18 


44.25 


8.16 


4 


S07 


23 09 46.14 


+00 00 48.9 


0.3518 


45.47 


2.10 


0.45 


44.80 


8.53 


2 


S08 


23 59 53.44 


-09 36 55.5 


0.3585 


45.23 


1.22 


0.40 


44.55 


8.10 


4 


S09 


00 59 16.11 


+15 38 16.1 


0.3542 


45.37 


0.71 


0.22 


44.78 


8.13 


4 


SIO 


01 01 12.07 


-09 45 00.8 


0.3506 


45.48 


1.02 


0.27 


44.93 


8.25 


4 


Sll 


01 07 15.97 


-08 34 29.4 


0.3557 


45.34 


0.52 


0.14 


44.85 


8.00 


4 


S12 


02 13 40.60 


+13 47 43.3 


0.3570 


45.43 


0.97 


0.28 


44.62 


8.67 


4 


S21 


11 05 56.18 


+03 12 43.3 


0.3534 


45.88 


2.15 


0.34 


44.06 


8.79 


1 


S23 


14 00 16.66 


-01 08 22.2 


0.3510 


45.42 


1.11 


0.29 


44.75 


8.70 


3 


S26 


15 29 22.26 


+59 28 54.6 


0.3691 


45.13 


0.52 


0.27 


44.38 


8.02 


4 


S27 


15 36 51.28 


+54 14 42.7 


0.3667 


45.23 


0.95 


0.36 


44.74 


8.10 


4 


S28 


16 11 56.30 


+45 16 11.0 


0.3680 


45.13 


0.76 


0.03 


44.66 


7.90 


4 


S29 


21 58 41.93 


-01 15 00.3 


0.3575 


45.43 


0.59 


0.10 


44.54 


7.94 


4 


S31 


10 15 27.26 


+62 59 11.5 


0.3504 


45.38 


0.29 


0.08 


44.58 


7.94 


2 


SSI 


08 04 27.99 


+52 23 06.2 


0.3566 


45.41 


0.39 


0.11 


44.89 


7.75 


4 


SS2 


09 34 55.60 


+05 14 09.1 


0.3675 


45.21 


0.33 


0.13 


44.54 


7.72 


3 


SS9 


12 58 38.71 


+45 55 15.5 


0.3704 


44.92" 


0.93 


0.26 




8.05 


4 


SSIO 


13 34 14.84 


+1142 21.5 


0.3658 


45.30" 


2.26 


0.43 




7.94 


4 


SSll 


13 52 26.90 


+39 24 26.8 


0.3732 


44.66" 


0.51 


0.15 




8.11 


4 


SS12 


15 01 16.83 


+53 31 02.4 


0.3626 


45.46" 


3.24 


0.52 




8.15 


4 


SS13 


15 05 41.79 


+49 35 20.0 


0.3745 


44.94" 


0.98 


0.38 




7.63 


4 



NOTE. — Col. (1): Object name. Col. (2): RA. Col. (3): DEC . Col. (4): redshift Col. (5): AGN bolometric lumino.sity measured 
from multi-wavelength data, adopted from Szathma rv et alj 120 111). Col. (6): AGN monochromatic luminosity at 5100A in 10"^^ erg 
s"^ units, corrected for host galaxy starlight based on the HST images, taken from Treu et al. i 2007); Benneit et al. 1 2010). Col. (7): 
AGN traction in t he observe d flux at 5100A, adopted from Treu et al. ( 2007); Bennert et al. 1 2010). Col. (8); AG N X-ray luminosit y 
from'Szathmai-v etaTI t20TTh . Col. (9); Black hole mass estimated from H/3 line width and L5100 adopted from lWoo et all 1200(1) ; 
[Bennert et al. 1 201^Col. (10); Number of pointings for AKARI observation 

^ AGN bolometric luminosity is calculated from optical luminosity (L5100) assuming that L/,^,/ = 9 X A L5100. 



PAH emission lines in various sources and found no system- 
atic difference in the measured PAH fluxes while Smith et alj 
(|2007) applied Drude profile to mid-IR spectra obtained from 
the Spitzer IRS. We chose a Gaussian profile, which is gener- 
ally adopted for the 3.3 /im PAH emission lines observed in 
starburst galaxies (Type-A sources in Tokunaga et al. 1991). 
In practice, we applied a Gaussian profile with a fixed width in 
the rest-frame for modeling the 3.3 /im PAH line. We adopted 
21 nm as a fixed width of the line profile, which has been de- 
termined from the PAH line profile of M82 (Tokunaga et al. 
1991) and used for our previous studies of the 3.3/im PAH 
line based on the AKARI IRC spectra (Imanishi et al. 2008, 
2010). 

Continuum subtraction is the largest source of uncertainty 
in measuring PAH emission line flux since absorption fea- 
tures are abundant in the near- and mid-IR spectral ranges. 
The vicinity of the 3.3 jim PAH feature is not an exception: 
There is a plethora of absorption features in the rest-frame 2 
- 5 pm range. We applied a linear fit to the continuum in 
the 3.0 - 3.6 /im range with two-pixel binning, after masking 
the location of the PAH feature around the 3.3 /im, as simi- 
larly practiced in other studie s (e.g. Imanishi & Wada 2004; 
llmanishi et al]|2010l:IOi et"an i2010: Imanishi et al. 2011). 

While fitting the 3.3 /im PAH feature by a Gaussian profile, 
we determined the peak and the central wavelength of the pro- 
file based on the minimization. Then, the line flux and the 
equivalent width were determined from the fit (see Table |2]l. 
The typical error of the measured flux density is ^ 3 x 10"'"* 



erg s"' cm"^ /tm"' . 

We detected the 3.3 /im PAH feature from 7 out of 26 ob- 
served galaxies, resulting in a detection rate of 27%. Com- 
pared to higher luminosity AGN, the 3.3 /im PAH detection 
rate of our lower-luminosity AGN is slightly higher although 
the difference is not easy to interpret owing t o the different 
exposu re time and flux limits. For example, llmanishi et alj 
( 1201 II) investigated the PAH emission features of 30 PG 
QSOs based on the ground-based observations, and detected 
the 3.3/im PAH emission from 5 QSOs, resulting in a ^ 
17% detection rate. If we exclude S3 and SS12, for which 
3.3/im PAH emission is marginally detected, then the detec- 
tion rate of our lower-luminosity AGN decreases to 19%, sug- 
gesting that the detection rate is similar to that of PG QSOs. 

Compared to the detection rate of other PAH lines at longer 
wavelength (e.g., 6.7, 7.7, 8.6, 11.2 /im) in local Seyfert 
galaxies, the detection rate of the 3.3/im PAH line in our sam- 
ple seems to be much lower (see Sales et al. 2010). This 
difference may be due to the fact that 3.3/im PAH is much 
weaker than other lines at longer wavelength, thus more dif- 
ficult to be detected. For better understanding of the relation 
among various PAH lines, spectroscopic studies with a large 
spectral range including 3.3/im and other PAH lines are re- 
quired. 

The detected 3.3 /im PAH line luminosity (L3.3 pah) ranges 
from ^ 0.9-5.7 x 10"*^ erg s~' while the upper limit of the 
PAH luminosity of the undetected galaxies is close to ~ 10"^^ 
erg s"'. Compared to the 3.3 /im PAH luminosity range of 
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Fig. 1 . — AKARI spectra of moderate-luminosity AGN at z~0.4 with the 
3.3 fim PAH emission detection. The location of the 3.3 ^tm PAH emission 
is indicated by an aiTow in each panel. 



Type I AGN host galaxies measured from ground-based ob- 
servations (L3,3 PAH=10^'' - 10"^^ erg s"'), the 3.3 ^m PAH lu- 
minosities from the AKARI slit-less spectroscopy are more 
than a factor of 10 higher. This is probably due to the aperture 
effect, reflecting the difference between nuclear and global 
star formation. Ground-based observations generally used a 
narrow-slit to probe nuclear starburst while AKARI slit-less 
spectroscopy provides the global star formation rate over en- 
tire galaxies. A direct comparison between ground-based and 
AKARI observations supports this interpretation. For exam- 
ple, a recent AKARI study by Imanishi et al. (2010) reported 
the 3.3 fim PAH luminosity of NGC 7469 as 6.3 x 10"^' erg 
s"' whil e the ground-based spectroscopy with a 1.6" slit 
width by llmani shi & Wada (2004) presented an upper limit 
of 2.7 X 10^°erg s"' , indicating that the nuclear PAH emission 
is a small fraction of the global PAH emission. 

For galaxies with no 3.3 PAH emission detection, we 
measured the upper limits. First, we averaged the 1 a errors 
at the wavelength range of the expected 3.3 fim PAH line. 
Then, we assumed a Gaussian profile with a peak flux, which 
is 3 times larger than the average 1 a error. In the case of the 
line width of the profile, we used a constant width of 21 nm 
in the rest-frame. The measurements and upper limits of the 
3.3 /xm PAH luminosity are listed in Table 2. 

3.2. Combined Spectra 

In order to construct the representative spectra of the sam- 
ple and to recover the undetected 3.3 /im PAH in individual 
galaxies, we constructed two combined spectra respectively 
for the entire sample and the galaxies with no 3.3 fim PAH 
emission detection. To construct stacked spectra, we used a 



composite non — detec:iori 

PAH 3,3 



Re'it wavelength (/lirn) 

Fig. 2. — Stacked spectra of the AKARI sample. The upper panel shows 
the spectrum stacked with the entire sample of 26 galaxies, while the bottom 
panel shows the spectrum spectrum stacked with the galaxies without the 3.3 
fim PAH emission detected. These stacked spectra are the sum of normalized 
individual spectra within the wavelength range of 3.1 to 3.5 fim in the rest 
frame. We also apply two pixel binning to the data points outside the 3.3 
fim PAH emission feature. 



wavelength range between 3. 1 - 3.5 /im in the rest frame and 
normalized individual spectra by a continuum fit. Then, we 
summed individual spectra without weighting. We also ap- 
plied two-pixel binning to the data points outside of the 3.3 
fim PAH line to smooth the continuum. The stacked spectra 
of the entire sample and the galaxies with no 3.3 fim PAH 
detection are presented in Figure |2] 

The combined spectrum of the entire sample clearly shows 
the 3.3 fim PAH emission. By fitting the PAH emission using 
a Gaussian profile with a width of 21 nm, we measured the 
EW of the 3.3 fim PAH fine. The EW of PAH in our sample 
is larger than that of QSOs. For ex ample, based on a sample 
of PG QSOs, llmanishi et aH (120111) reported that the EWs of 
the 3.3 fim PAH emission of individual galaxies are typically 
less than 10 nm while the EW of the 3.3 fim PAH emission 
measured from their combined spectrum is only 1 nm. The 
smaller EW of PAH emission in higher luminosity AGN may 
be interpreted to mean that PAH lines are diluted by higher 
AGN continuum. However, a larger sample is required to re- 
veal the origin of the trend. 

4. STARBURST-AGN CONNECTION 

In this section, we investigate how starburst and AGN activ- 
ity is related. First, we compare AGN activity with the global 
star formation rate measured from the AKARI observations 
(§ 4.1). Then, we investigate the correlation between nuclear 
starburst and AGN activity based on the collected data from 
the hterature (§ 4.2). 

4.1. Global PAH emission and AGN activity 

To investigate the relation between global star formation 
and AGN activity, we compare the luminosity of the 3.3 
fim PAH emission (L3 3 pah) with various AGN properties, 
i.e., optical. X-ray, and bolometric luminosities. First, in Fig- 
ure 3 (left) we compare the L3.3 pah with AGN continuum 
luminosity at 5100A (L5100) measured from high S/N optical 
spectra obtained with the Low Resolution Imagi ng Spectro- 
graph at the Keck Telescope tWoo et ahi2006; .Bennert et alj 
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TABLE 2 

Global 3.3 fiM PAH emission measurements based 

ON THE AKARI OBSERVATIONS 



KJ Uj CL L 


13. 3 


L3.3^m 


^*V3,3^m 




(lU erg s cm ) 


(10 erg s 


') (nm) 


(1) 


(2) 


(3) 


(4) 


SOI 


< 5.09 


<2.19 




S02 


< 10.47 


<4.36 




803 


3.10 


1.33 


9 


S04 


< 1.55 


<0.66 




SOS 


< 2.00 


<0.83 




S06 


< 2.03 


<0.93 




S07 


< 2.76 


<1.13 




SOS 


3.80 


1.63 


87 


S09 


< 31.41 


< 13.07 




SIC 








Sll 


4.40 


1.85 


21 


S12 


< 1.77 


< 0.75 




S21 


< 12.73 


<5.27 




S23 


< 1.97 


<0.80 




S26 


12.50 


5.73 


206 


S27 


<1.69 


<0.76 




S28 


6.65 


3.03 


43 


S29 


< 1.51 


<0.64 




S31 


< 1.43 


<0.58 




SSI 


< 1.77 


<0.74 




SS2 


< 1.84 


<0.83 




SS3 


< 1.57 


<0.66 




SS9 


< 1.72 


<0.79 




SSIO 


< 1.56 


<0.70 




SSll 


7.15 


3.36 


77 


SS12 


2.00 


0.88 


3 


SS13 


< 1.74 


<0.82 




Stacked 






21 


Note. — 


Col. (I): Object name. 


Col. (2): 3.3 /^m PAH line flux. 



Col. (3). 3.3 fj,m PAH line luminosity. Col. (4). Equivalent width of the 
3.3 fim PAH emission line. 



120101) . Although the optical luminosity range is rather limited, 
we find no strong positive correlation. Instead, a weak neg- 
ative correlation between L3 3 pah and the optical luminosity 
seems to be present among PAH-detected objects. Assuming 
that the strengths of the 3.3 and 7.7/im PAH emissions are cor- 
related, this result is inconsistent w ith previous studies (e.g., 
iNetzer et al.l2007l:lLutz et al.120081) . which reported a positive 
correlation between 7.7/xm PAH luminosity and 5100A lumi- 
nosity. One potential bias is the contribution from host galaxy 
starlight to L5100 since for low luminosity AGN, such as the 
sample considered here, the amount of starlight c an be com- 
parable to t he AGN flux in the optical range (e.g.. lWoo et al.l 
|2()06; Benn ertet al.] |2010). 

To correct for the host galaxy contribution, we mea- 
sured the nuclear luminosity at 5100A based on the 2- 
dimensional AGN-host galaxy decomposition analysis using 
high-resolution images obtained with the Hubble Space Tele- 
scope (for details, see Treu et al. 2007; Bennert et al. 2010). 
The mean AGN-to-total flux ratio of our sample is 0.26 (see 
Table 1), indicating that the correction for the host galaxy con- 
tribution at 5100A is necessary. After correcting for the host 
galaxy contribution, we compare the pure AGN optical lumi- 
nosity at 5100A with the PAH luminosity as shown in Figure 
3 (right). The nuclear luminosity (Lsioo.nuc) becomes signifi- 
cantly smaller than the observed luminosity at 5100A. How- 
ever, this correction does not significantly change the trend 
between L3 3 pah and AGN continuum luminosity. 

These results are rather contradictory to the previous stud- 




44 45 43 



Log L5100A (erg s-i) Log L^iqo^^^^ (erg s"') 



Fig. 3. — Comparison of the 3.3/im PAH luminosity (L3 3 pah) with the 
observed total luminosity at the rest-frame 5100 (L5100) (left) and AGN lu- 
minosity corrected for host galaxy starlight contribution (right). Filled circles 
represent the PAH detected galaxies while open circles represent upper lim- 
its of L3 3 PAH. Note that the monochromatic luminosity at 5100A includes 
the contribution from starlight, which is substantial fraction of the total flux 
(74% on average). 

ies that claimed positive correlations between AGN activity 
and various s tarburst indicators, i.e., L3 3 pah and far-IR lu- 
minosity (e.g.'N etzer eTaPIIOOTHWatabe et allllOOllOi et alJ 
I2OIO). Note that there are many Type 1 AGN witiiout PAH 
detections in the previous studies as well as our observation. 
By including these upper limits, the correlation between PAH 
and AGN luminosities can substantially weaken. On the other 
hand, since the luminosity range of our sample is rather lim- 
ited, we cannot conclude whether or not there is a correlation 
between L3 3 pah and Lsioo in a larger luminosity range. 

Instead of optical luminosities, we can use X-ray luminos- 
ity, which is a good indicator of AGN activity since X-ray 
is dominantly radiated from accretion disk in Type I AGN. 
In Figure 4 (left), we compare L3.3 pah with AGN X-ray lu- 
minosity integrated over the 0.5 to 8 keV range, measured 
from Chandra X-ray images (Szathmary et al. 2011). Within 
the order of magnitude range of the X-ray luminosity, we see 
no strong positive correlation. We also used the bolometric 
luminosity of AGN determine d from multi-wavele ngth SED 
fitting analysis (for details, see Szathmary et alj201 li) to com- 
pare with L3 3 PAH- As shown in Figure 4 (right), we see no 
strong correlation between L3 3 pah and the AGN bolometric 
luminosity within the hmited range of the observed luminos- 
ity. 

4.2. Nuclear PAH emission and AGN activity 

In this section, we investigate how the nuclear 3.3/im PAH 
emission is related to AGN activity. We collected from 
the literature all 3.3/im PAH measurements of Type I 
AGN, for which AGN luminosity and black hole mass 
can be compared. We found 3.3/im PAH measure- 
ments for 20 local Seyfert 1 galaxies and PG QSOs 
based on the previous ground-based long-slit observations 



6 



TABLE 3 

Nuclear 3.3/iM PAH luminosity and AGN properties 



Object 


z 


log LpAH3.3 
(erg/s) 


Ref. 


log Mbh 

(M©) 


Method 


Ref. 


log AL5ioo,n„c 

(erg/s) 


Method 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



NGC 3227 


0.0039 


39.36 


2 


7.60 


reverberation 


2 


42.48 


bulge-disk decomposition 


2 


NGC 4235 


0.0008 


39.64 




7.71 


single-epoch virial 


this work 


42.08 


from Lhq 


this work 


NGC 4748 (MCG-2-33-34) 


0.0146 


40.18 




6.39 


reverberation 


1 


42.74 


from spectral fitting 


1 


NGC 5273 


0.0036 


39.30 




6.85 


single-epoch virial 


this work 


41.45 


from Lhq 


this work 


NGC 5940 


0.0339 


40.80 




7.95 


single-epoch virial 


this work 


43.37 


from Lhq 


this work 


Mrk 335 


0.0258 


40.49 




7.13 


reverberation 


2 


43.78 


bulge-disk decompositoin 


2 


Mrk 509 


0.0344 


41.00 




8.14 


reverberation 


2 


44.16 


bulge-disk decompositoin 


2 


Mrk 530 


0.0295 


40.88 




8.06 


Mbh"""* relation 


3 


43.71 


from Lboi. 


3 


Mrk 618 


0.0356 


40.85 




8.34 


MBH-Lbulgc relation 


4 


44.03 


from Mbh and Eq. 1 


4 


Mrk 766 (NGC 4253) 


0.0129 


40.32 




6.23 


reverberation 


1 


42.81 


from spectral fitting 


1 


Mrk 817 


0.0315 


41.00 




7.67 


reverberation 


2 


43.64 


bulge-disk decomposition 


2 


3C 120 


0.0330 


41.04 




7.72 


reverberation 


2 


44.09 


bulge-disk decomposition 


2 


PG 0157+001 (Mrk 1014) 


0.1631 


42.32 


3 


8.31 


single-epoch virial 


5 


44.89 


spectral fitting 


5 


PG 1211+143 


0.0809 


41.43 


3 


8.14 


reverberation 


2 


44.70 


bulge-disk decomposition 


2 


PG 141 1+442 


0.0896 


41.49 


3 


8.63 


reverberation 


2 


44.52 


bulge-disk decomposition 


2 


PG 1416-129 


0.1289 


41.40 


3 


8.83 


single-epoch virial 


5 


44.55 


spectral fitting 


5 


PG 1440+356 


0.0791 


41.51 


3 


7.80 


single-epoch virial 


6 


44.54 


spectral fitting 
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NOTE . — Col. (I): Obj ect name Col. ( 2): Redshift. Col. (3): 3.3 /jm PAH luminosity Col. (4): reference for L3.3 pah- 1 . Ilmanishi & Wad j tlOOl . 2. iRodriguez-Ardila & Viegaj 120031) 
3.'Imanishi et al. ( 201 IJ). 4. IOi et alj 12010) . 5. Ilmanishi et alj 12010) . Col. (5): Black hole mass Col. (6): reference for black hole mass I. Bentz et al. (200%). 2. Peterson et al. (2004). 3. 
Nelson & Whittle (1995). 4. Ryan et al. (2007). 3. Col. (7): AGN monochromatic luminosity atSIOOA, Col. (8): reference for L5100. 1. Park et al. (201 1). 2. Bentz et al. (2009a). 3. Woo & 
Uny (2002). 4. Ryan et al. (2007). 5. Ho & Kim (2009). 6. Shang et al. (2007). 
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Fig. 4. — The correlation of Lj jpah with X-ray luminosity (io,5-8A:ev) mea- 
sured from Chandra X-ray images. Symbols are same as in Figure[3] 



Fig. 5. — The correlation of nuclear L3 3 pah with AGN optical luminos- 
ity. Filled squares represent L3 3 pah measurements from ground-based spec- 
troscopy with a narrow-slit while open circles represent the estimated nuclear 
L3.3 PAH of the AKARI sample, after dividing by a factor of ~23, which is 
the ratio between global and nuclear' PAH luminosity of NGC 7469. 



jRodriguez- Ardila & Viegas 2003 ; Ilmanishi & Wadal l2004t 
IWatabe et al.i i2008i: lOi et al.i i20l'd Ilmanishi et alJ l20nl) . 
These measurements based on the long-sht spectroscopy with 
a narrow slit represent the strength of the nuclear 3.3/im PAH 
emission. Except for 3 AGN without necessary measurements 
of AGN properties in the literature, we determined black hole 
mass and AGN continuum luminosity at 5100A. Thus, we 
have a sample of 17 local Type I AGN as listed in Table 3. 

For black hole masses, we adopted the reverberation map- 
ping measurements from Peterson et al. (2004) and Bentz 



et al. (2009b), after calibrating with the mo s t rece nt deter- 
mination of the virial factor from ' Woo et a D dMol) . When 
reverberation mapping results were not available, we used 
the single-epoch virial method as used for the moderate- 
luminosity AGN sample. In practice we use the following 
equation to estimate black hole masses (Park et al. 201 1): 



Mbh = lO'-^^'Mo 



V 1000 km s"' 



XL 



6l00jmc 



1044 erg s- 



0.518 



(1) 
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where a^p is the line dispersion of the H/3 line and Lsioo.nue 
is the AGN continuum luminosity at 5100A corrected for the 
host galaxy starlight contribution. For NGC 4235, NGC 5273, 
and NGC5940 with no line width measurements in the liter- 
ature, we fit the broad Ha emission line and measured the 
line width and luminosity using the spectra from the Sloan 
Digital Sky Survey. For two objects (Mrk 530 and Mrk 618) 
without any available spectra, we used the black hole mass 
-galaxy propert y relations, i .e., Mbh- errand MfiH-Lbuige re- 
lations (Gtiltek in et al.ll2009h to calculate black hole masses 
(see Table 3 for details). 

In the case of AGN power, we collect AGN continuum lu- 
minosity at 5 100 A as a proxy for the bolometric luminosity of 
AGN. For high luminosity QSOs, i.e., PG QSOs, it is not nec- 
essary to correct for the host galaxy starlight contribution. In 
contrast, for low luminosity AGN, Lsioo should be corrected 
for the host galaxy contamination. In general we obtained 
Lgioo.nuc measurements from previous studies, which were 
based on either bulge-disk decomposition analysis using HST 
images or spectral fitting analysis using high quality spectra. 
For 3 objects, namely NGC 4235, NGC5273, and NGC 5940, 
we determined AGN continuum luminosity from Ha line lu- 
minosity using the tight correlation between Lsioo,mic and Lua 
(see Eq. 1 and 2 in Greene & Ho 2005). 

In contrast to the ground-based observations, L3 3 pah mea- 
sured from AKARI observations represents the global star- 
burst as discussed in § 3.1. Thus, in order to include the 
PAH emission measurements from AKARI observations, it 
was necessary to calibrate them for representing nuclear star- 
burst. The relation between nuclear and global starburst for 
Type I AGN has not been studied in detail (Imanishi et al. 
I2OIOI) . Although the nuclear-to-global starburst ratio can vary 
for individual objects, we used the case of NGC 7469 to es- 
timate the difference. As discussed in § 3.1, NGC 7469 was 
observed with a ground-based spectrograph with a narrow-slit 
and the AKARI IRC, respectively. The flux ratio between nu- 
clear L3 3 PAH and the upper limit of the global L3 3 pah is 
^0.04. Thus we took this ratio as a calibration factor and 
scale down the AKARI measurements by a factor of 0.04. 

In Figure 5, we present the correlation between the nu- 
clear 3.3/Ltm PAH luminosity and AGN optical luminos- 
ity at 5100A for the combined sample. The Type I 
AGN from ground-based observations clearly show a corre- 
lation, suggesting the nuclear starburst- AGN connection. The 
L3.3 PAH of the moderate-luminosity AGN at z~0.4, which are 
calibrated for representing the nuclear PAH emission, are con- 
sistent with the trend found by nuclear L3 3 pah measurements 
although they are more scattered probably due to the uncer- 
tainty of the calibration. 

Using the nuclear L3 3 pah measurements from ground- 
based spectroscopy, we fit the relation between L3 3 pah and 
AGN luminosity. Since most of the L3 3 pah measurements in 
the literature do not have error estimates, we fit the best-slope 
without accounting for the measurement uncertainties. The 
best-fit relation is: 



log( 



-PAH3.3 



1041 



ergs- 



■j) = 0.06 ±0.01 +(0.77 ±0.01) xlog(-^^^i^). 

10^^ ergs ' 

(2) 



The best-fit slope of ~0.8 indicates that the correlation is not 
linear This implies that for higher luminosity AGN, the nu- 
clear starburst-to-AGN luminosity ratio is smaller, albeit with 
the small sample size. 
In Figure 6, we compare nuclear L3.3 pah with the Edding- 
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Fig. 6. — Correlation of nuclear' L3 3 pah with AGN optical luminosity. 
Filled squares represent L3 3 pah measurements from ground-based spec- 
troscopy with a narrow-slit while open circles represent the estimated nuclear 
L3.3 PAH of AKARI sample, after dividing by a factor of ~23, which is the 
ratio between global and nuclear PAH luminosity of NGC 7469. 



ton ratio and black hole mass to test whether the accretion 
rate and black hole mass correlates with starburst activity. We 
estimate the Eddington ratio by multiplying a factor of 9 to 
the AGN optical luminosity. In the case of the Eddington ra- 
tio, there is a general trend that higher Eddington ratio objects 
show higher nuclear starburst activity although there is a large 
scatter Black hole mass shows a similar trend with L3 3 pah in 
that the host galaxies of higher mass AGN show higher nu- 
clear starburst luminosity. This can be interpreted as the trend 
that higher mass AGN in the sample are generally higher lu- 
minosity AGN, i.e., PG QSOs. Thus, the relation between 
black hole mass and PAH luminosity reflects the correlation 
between AGN power and L3 3 pah- 

5. DISCUSSION AND CONCLUSIONS 

To investigate the starburst- AGN connection, we compared 
the global and nuclear 3.3 /im PAH emission with AGN lu- 
minosity using a sample of Type I AGN. The global 3.3 
/im PAH emission luminosities are measured for moderate- 
luminosity AGN at z^O.4, based on our new AKARI slit-less 
spectroscopy, which covers the entire host galaxies. We find 
no strong correlation between 3.3 /im PAH and AGN optical 
luminosities within the limited AGN luminosity range lO^^ - 
10^^ erg s"'. This result may imply that the global star for- 
mation is not strongly related with AGN activity although we 
cannot rule out a weak correlation with a considerable scatter. 
We note that since the AGN luminosity range is very limited, 
whether the global star formation rate correlates with AGN 
activity is inconclusive. 

In contrast, when we compare AGN luminosity of the lo- 
cal Seyfert 1 galaxies and PG QSOs with their nuclear 3.3 
lim PAH emission luminosity, we detect a strong correlation, 
suggesting that that AGN activity is related to the nuclear star- 
burst. These results are consistent with the findings of previ- 



8 



ous studies on the starbur st-AGN con nection ('Imanishi"2003'; 
Imanishi & Wada 2004; Netzer etalj 1200 7: Lutz et al. 2008; 
Oi et alJl201Ql) . 

The slope of the L3 3 pah correlation with AGN luminos- 
ity in Figure 5 shows that the L3.3 pah to AGN luminosity ra- 
tio decreases toward higher AGN luminosity. At face value, 
it implies that the nuclear starburst activity is slightly sup- 
pressed in the host galaxies of high luminosity AGN (for dif- 
ferent predictions on the starburst-to-accretion ratios, see e.g., 
Kawakatu & Wada 2008 ; Ballantyne 2008). A similar trend 
was noted by iLutz et all (i2Q08) in comparing 7.7 fim PAH 
emission and AGN luminosities. In their study, high lumi- 
nosity AGN with L5100 > 10"*^ erg s~^ show a much l ower 
starburst-to-AGN ratio (see Figure 6 in lLutz et al.ll2008l ). im- 
plying that the black hole growth rate normalized by the star- 
burst rate is a factor of 10 higher in their high redshift QSOs 
than in the local AGN. In our local sample the decrease of 
the starburst-to-AGN ratio with increasing AGN luminosity 
is not as dramatic as in high redshift QSOs, but if confirmed. 



it implies that for low-luminosity AGN the growth of galaxy 
centers is faster than black hole growth in the present-day uni- 
verse. However, it is not clear whether black hole growth is 
slower or faster than the total host galaxy growth since only 
nuclear starbursts are compared in our analysis. 

We note that there are many AGN with no PAH emission 
detection. Including the upper limits of these objects signifi- 
cantly weakens the correlation between L3.3 pah and AGN lu- 
minosity. Although the true PAH luminosities of these objects 
may be close to the upper limits within a factor of few, a de- 
tailed investigation using a large sample with lower flux limits 
is necessary to probe the slope and scatter of the correlation. 

This work has been supported by the Basic Science Re- 
search Program through the National Research Foundation 
of Korea (NRF) funded by the Ministry of Education, Sci- 
ence and Technology (MEST), No. 2010-0021558. M.I. ac- 
knowledges support by Grants-in-Aid for Scientific Research 
no. 22012006. This research is based on observations with 
AKARI, a JAXA project with the participation of ESA. 
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